Solar supergranulation 
as a wavelike phenomenon 



<sj ■ L. Gizon*, T. L. Duvall Jr.t, and J. Schou* 

O ■ 

* W. W. Hansen Experimental Physics Laboratory, Stanford University, Stanford, CA 94305, USA 

' Laboratory for Astronomy and Solar Physics, NASA Goddard Space Flight Center, Greenbelt, MD 20771, USA 

ri ' Supergranulation on the surface of the Sun is an organized cellular flow pattern with 
a characteristic scale of 30 Mm[[l]]. It is superficially similar to the well understood 
^ . granulation [§] that operates at the 1.5 Mm natural scale of convection, which has 
led to the conventional view that supergranulation has its origin in the convective 
motion of cells of gas[^, ||, though this does not explain the observation that the 
supergranulation pattern appears to move faster around the Sun than the bulk of 
its surface^], §], |7|- A wave origin has been proposed[]8|] for supergranulation that 
qq . may explain the superrotation, but it has never had much support. Here we report 
O ! that the supergranulation pattern has oscillatory components with periods of 5-10 
^ \ days, for which the best explanation is a spectrum of traveling waves. We show that 
' there is excess power in the prograde and equatorward directions, which explains the 
Qlt observation of superrotation. 

To study supergranulation, we use a 60-day sequence of line-of-sight Doppler velocity images of the 
^ \ Sun's surface, obtained by the Michelson Doppler Imager on board the Solar and Heliospheric 
Observatory (SOHO). We apply the technique of time-distance helioseismology |T0[| to obtain maps 



of the horizontal divergence of the flows just below the solar surface. Unlike raw Doppler images, 
^ the divergence signal has uniform sensitivity across the solar disk and has few systematic problems. 
Every 6 hr, an 8 hr section of data is used to obtain a divergence map. This is done by measuring the 
time it takes for surface gravity waves (f modes) to propagate horizontally between a given point 
on the solar surface and a concentric annulus with radius 8 Mm. The difference in travel times 
between inward and outward propagating waves is a proxy for the horizontal divergence of the flow 



field in the first 2 Mm beneath the surface. The details of this analysis are described elsewhere |TT 
Supergranules appear as cellular patterns of horizontal outward flow in the divergence maps. 

We considered five latitude bands, 90° in longitude by 10° in latitude, centered at latitudes 
0°, ±15°, and ±30°. For each latitude band, the divergence signal is Fourier transformed in three 
dimensions to make power spectra as a function of angular frequency, u, and horizontal wavevector, 
k. Figure 1 shows a cylindrical cut in a power spectrum at a constant wavenumber, k = ||k||, typical 
of the supergranulation. The power peaks at a non zero frequency implying that the supergranular 
divergence field oscillates from positive to negative values. This observation alone argues against a 



simple model of convection JT2], [Tj|. Figure 1 is characteristic of waves propagating in all directions 
at a well defined frequency. The variation of power with the direction of propagation, <ft, indicates 
a strong anisotropy. We thus suggest that we are observing a superposition of traveling waves, and 
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show below that this hypothesis provides a solution to the superrotation puzzle. 

Waves that propagate horizontally through a moving fluid have a dispersion relation of the form 
uj = uj (k) + k • v, where Uq is the wave frequency in the co-moving frame of reference and v 
is the advection velocity. At fixed wavenumber, we thus interpret the sinusoidal variation of the 
frequency with <fi to be due to a horizontal flow (Fig. 1). This technique for measuring material flows 
is equivalent to the ring-diagram technique [111] used in helioseismology. Figure 2 shows the inferred 



rotation and meridional circulation, both similar to that of the small magnetic features ||15|, jig] 



The rotation corresponds to the helioseismic value at a depth of 20 Mm, suggesting that the waves 
probe the region of increased rotation just beneath the solar surface. Correlation tracking of the 
supergranulation pattern, on the other hand, gives a larger rotation rate, larger that the rotation 
of the plasma measured at any depth in the solar interior 0, as well as a puzzling equatorward 
meridional motion (Fig. 2). 

The anomalous results from correlation tracking can be understood as a consequence of the 
pronounced anisotropy of wave power (Fig. 3). The flow velocities inferred from the correlation 
tracking are effectively a power-weighted average of the true advection velocity, v, and the wave 
propagation speed (about 65 m s _1 ). The supergranulation pattern appears to superrotate because 
wave power in the prograde direction is almost twice as large as in the retrograde direction. Since 
the change in the rotation over the outer 50 Mm is a significant fraction of the phase speed, a 
buried source of excitation could introduce an anisotropy. Note also that a model of standing-wave 
oscillatory convection advected by a flow would fail to reproduce the anisotropy in wave power. 

We have shown that the power spectrum of the divergence field of the supergranultion is con- 
sistent with an anisotropic spectrum of traveling waves. One possibility |J for the nature of the 
waves is that they are driven by the Coriolis force, analogous to the oceanic Rossby waves fl8H and 
perhaps related to the recently observed solar r modes|13|. However, wave periods are found to be 
significantly smaller than the rotation period (Fig. 4) and we found no evidence for a latitudinal 



variation in the dispersion relation. Another possibility is traveling-wave oscillatory convection) 19 



For example, numerical simulations of turbulent convection in oblique magnetic fields p(J exhibit 
solutions that take the form of traveling waves, where the tilt of the convection cells, their wave 
speed, and direction depend on the field strength and obliquity. Finally, we note that the anisotropy 
in wave power is important as it could provide a mechanism of angular momentum transport in the 



upper convection zone, that could help explain the variation of rotation with latitude[21 
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Figure 1 Cylindrical cut at constant wavenumber in the 3D power spectrum of the divergence 
signal near the equator, with k = 120/ R where R is the solar radius. Black indicates larger power. 
The azimuth, 0, is the angle between k and the direction of solar rotation (West). The top ridge 
is the power in the waves traveling in the direction given by <fi. The bottom ridge is for waves 
traveling in the opposite direction and is identical to the upper ridge, but shifted by 180°. We 
applied an MTF correction estimated from the data. To measure the flows, we first fit, at each <j>, 
the sum of two Gaussian functions in frequency with independent central frequencies and common 
line widths. We then fit curves (double lines) of the form u = ±uo + kv x cos <j) + kv y sin to the 
central frequencies, where ojq is the wave frequency, v x is the velocity in the direction of rotation, 
and v y is the meridional circulation. For each latitude, v x and v y are measured to be roughly 
independent of k in the range 50 < kR < 200, consistent with the interpretation as a flow. The 
raw data was tracked at the Carrington rotation rate, leading to a small prograde frequency shift 
for the equatorial band. 
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Figure 2 Flows inferred from the advection of the supergranulation waves, a, Solar rotation (di- 
amonds). The magnetic feature rate|15j (solid line) is very similar and the photospheric rate|| 
(dashed) is lower. The squares show the results using correlation tracking of the supergranula- 
tion pattern with a 24 hr delay, consistent with an earlier estimate || (dotted), and significantly 
higher than rotation inferred at any depth in the solar interior^], b, Northward meridional flow 
(diamonds). The magnetic feature ratcjl6j (solid) is again similar. For comparison, the dashed 
line shows a time-distance measurement |ll 7|| , while the squares show the anomalous results obtained 
using correlation tracking with a 24 hr delay. 
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Figure 3 Anisotropic wave power, a, Directional distribution of power at latitude —30° (solid 
line). The azimuth of maximum power (cross) is measured by fitting a sine (dashed), b, Azimuth 
of maximum power versus latitude, with 1-a error bars. Wave power is maximum in the direction 
of rotation and toward the equator in both hemispheres. The dashed line is drawn for the azimuth 
equal to the negative of the latitude. 
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Figure 4 Power spectrum corrected for rotation and meridional circulation and averaged over 
azimuth. Contour lines are drawn at 15%, 30%, 45%, and 75% of the maximum value. The 
diamonds show the frequencies, ivo(k), deduced from the fits (see Fig. 1). The empirical dispersion 
relation u 2 = ak with a = 7.5 x 10 -4 m s" 2 is overplotted for reference (thick line). We did not 
find a significant latitudinal variation in the dispersion relation, for latitudes less than 35°. The 
relatively large spread of the power indicates that the waves have a low quality factor. The width 
of the power increases rapidly with k, presumably due to the increased efficiency of the damping 
on small scales. Note that the distribution of power as a function of frequency is only affected by 
the known temporal window function, while the wavenumber dependence includes effects of the 
telescope optics and of the time-distance analysis. 
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